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ABSTRACT 

Reduced models are matrix integrals believed to be related to the large N limit of gauge 
theories. These integrals are known to simplify further when the number of matrices D 
(corresponding to the number of space-time dimensions in the gauge theory) becomes large. 
Even though this limit appears to be of little use for computing the standard rectangular 
Wilson loop (which always singles out two directions out of D), a meaningful large D 
limit can be defined for a randomized Wilson loop (in which all D directions contribute 
equally). In this article, a proof-of-concept demonstration of this approach is given for the 
simplest reduced model (the original Eguchi-Kawai model) and the simplest randomization 
of the Wilson loop (Brownian sum over random walks). The resulting averaged Wilson 
loop displays a scale behavior strongly reminiscent of the area law. 



1 Introduction 



Reduced models (for a textbook exposition, see [T], for a review of more recent develop- 
ments, see [2]) are matrix integrals whose large N limit is related to large N limits of gauge 
theories. The simplest and historically the first of these models, named after Eguchi and 
Kawai j3] , is simply obtained by a full dimensional reduction of a pure gauge theory to zero 
dimensions. The corresponding action reads 
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Here, the Greek indices label the D dimensions (of the original gauge theory space-time) 
and the Roman indices label the adjoint representation of the gauge group SU(N) with gen- 
erators t a . Various quantities are obtained by taking "path integrals" (which are ordinary 
integrals for this zero-dimensional case): 
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Of particular interest is the reduced model analog of the Wilson loop: 



ma*)) 



dA a „ — tr P exp 
M N 



L 



i / ds^{s)A a t a 



-S[A] 



Lo Lo 



dA«e~ s WJ2(i) n ds x ds 2 --- / ds B -tr ]J e k (s k )A^t a 



(4) 



Even though the Eguchi-Kawai model was originally hoped to reproduce large N prop- 
erties of the gauge theory from which it is obtained by dimensional reduction, it later turned 
out that this equivalence is undermined by a spontaneous breakdown of the R D symmetry 
that shifts the values of the matrices A by matrices proportional to unity. Variations of 
the Eguchi-Kawai model (the so-called twisted [1] and quenched [S] models) have been pro- 
posed to reproduce the large iV limit of the gauge theory (in particular, the reduced model 
Wilson loop of the type (jlj) would actually equal the gauge theory Wilson loop in the large 
N limit). However, they have been subsequently shown to suffer from non-perturbative 
problems [6]. A new generation of reduced models has been subsequently developed [21 E], 
and it is believed to reproduce correctly the properties infinite volume gauge theories. All 
these more sophisticated reduced models are closely related to the original Eguchi-Kawai 
model, even if computationally less tractable. 

In this paper, all derivations will be performed in the context of the Eguchi-Kawai model. 
The ideas are nevertheless quite general and should apply to more realistic reduced models 
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as well. This, however, would require a systematic development of the large D techniques 
(see below) for a broader class of reduced models. Likewise, introducing supersymmetry 
may protect the above-mentioned R D symmetry from spontaneous breaking, and would 
make the straightforward reduced model obtained by dimensional reduction of the gauge 
theory describe the large N limit (see [8], for recent considerations in the context of string 
theory, see [9]). 

A further significant simplification occurs in reduced models when the number of matri- 
ces D (corresponding to the number of space-time dimensions of the original gauge theory) 
becomes large. For the Eguchi-Kawai model, this simplification can be demonstrated by 
the following technique [TO]. (The large D techniques have been recently applied to matrix 
quantum mechanics [11] as well, though matrix integrals will remain our main focus in the 
present paper.) One first introduces an auxiliary field h ab , which is a real symmetric tensor, 
with the following action. 



S[A,h] = ^\ abcd (h ab h cd -A^h cd -h ab A c u A d u ) (5) 



= ^ abcd h ab h^ + ^K ab AlA b ^ (6) 
where K is a dimensionless kernel 

K ab = --\ abcd h cd (7) 
9 

for the Gaussian quadratic form of A a . Integrating out the auxiliary field, one reproduces 
the original action (CQ). 

Since (JBJ) is quadratic in A* it can be integrated out first. The propagator for A® is 
given by 



(A a A b v ) = gS^K^U. 



By rescaling the h as h ab = h ab /gvD, which is now dimensionless, we can rewrite the 
effective action as 

S cS = ^ {ttIuK + \\ abcd h ab h cd \ . (9) 



Therefore, in the large D limit, the integration over h is dominated by the saddle point and 
one can perform a systematic 1/D expansion. The saddle point equation can be given as 

(K-'U 6 -^ + X abcd ~h cd = 0. (10) 
oh ab 

We denote the relevant saddle point h^ b \ and if one assumes that it preserves the SU(iV) 
symmetry, h^ b = vd ab . From the saddle point equation (ITU1) . we obtain v = 1/V2N, i.e., 
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The saddle point values of K and its inverse (which is all one needs for computations to 
leading order in 1/D) are 



V2ND5 ab , 



1 



(12) 
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Then, from (JBJ), to compute the value of a certain observable to leading order in 1/D, one 
simply needs to contract all the A's with 



and sum over all such contractions. 

It is interesting to see what these techniques have to say about the expectation value 
of the Wilson loop (BJ, a crucial observable for gauge theories. It may come as a disap- 
pointment that there is no non-trivial large D limit for any fixed Wilson loop (i.e., Wilson 
loop for a given contour). A simple fixed Wilson loop singles out directions in the space 
of the A^s. On the contrary, in a theory with a large number of degrees of freedom (such 
as a large D limit of reduced models), simplifications occur for observables averaged over 
a significant fraction of these degrees of freedom, not for observables to which just a few 
degrees of freedom contribute. 

Mathematically, one can try to take a naive large D limit of the Wilson loop 01]). For 
simplicity, one can take the contour to be rectangular in the (1, 2)-plane. Then, only powers 
of A\ and A2 appear in the expansion of the path-ordered exponential in (jl]). According to 
( TT31 contractions of A\ and A2 will always produce factors of g/sJIND. On the other hand, 
organizing the large N and D limit as a power expansion ml/N and 1/D requires keeping 
the modified 't Hooft coupling g 2 ND fixed [10]. Hence, any combination of powers of A\ 
and A2 will be suppressed by powers of 1/D in the large D limit, and the D — > 00 limit of 
the simple rectangular Wilson loop will equal 1, containing no information on the dynamics 
of the theory. It is only when one considers combinations of the sort A^A^ that contracting 
the /i and v indices in ffl3]) lends an additional factor of D, and hence a non-trivial D — > 00 
limit. This is in agreement with the heuristic observations of the previous paragraph. 

There is a natural modification of the traditional Wilson loop observable that does not 
suffer from the lack of a meaningful and tractable large D limit. Namely, one can perform 
averaging of the Wilson loop over a large number of different contours. If the chosen 
ensemble of contours is rich enough, many directions (many different A^s) will contribute, 
and the large D limit will become dynamically non-trivial. In this paper, we shall consider 
the simplest ensemble of uncorrelated random walks (Brownian motion). However, the idea 
is quite general and should apply to more complex ensembles of loops as well. 

The Brownian randomized Wilson loops are certainly quite different from the regular 
Wilson loops familiar to the reader. For one thing, the typical fractal dimension [12] of 
a random walk trajectory is 2 (and it is nowhere differentiable). Still, one can give two 
reasons why considering sums over such ensembles of contours is a reasonable idea. First, (in 
a Euclidean theory) real trajectories of quantum particles (quarks) have fractal dimension 
2 (and not 1 like the idealized infinitely heavy quark implicit in the usual Wilson loop). 




(13) 
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Second, the difference between regular and random Wilson loops appears to diminish in 
a sense for a large number of dimensions (which is the case that concerns us here). To 
clarify this, one can consider a random walk with L steps on a D-dimensional hypercubic 
lattice. This walk would occupy a volume with LP I 2 points, whereas the walk itself contains 
L points, so (since there are no long range correlations) one can estimate the fraction of 
self-intersections (relatively to the total number of points) as L/L D ^ 2 . So, of course, in two 
dimensions self-intersections typically occupy a finite fraction of the entire trajectory and 
there is nothing like a notion of the perimeter and interior of the loop. Yet, in a very large 
number of dimensions, intersections become progressively more rare, and one should expect 
the random Wilson loops to look more and more tame. 

In what follows, we shall consider two cases of Wilson loop randomization. First, we 
shall generate the Wilson loop contour by a completely unconstrained random walk. In this 
case, the contour does not necessarily start and end at the same point, so a gauge theory 
interpretation of such a 'loop' is questionable (as it is not gauge invariant). However, the 
computation is particularly simple for this case, hence it is a nice toy model to develop the 
relevant techniques. Thereafter, we shall consider a sum over closed Wilson loops generated 
by a random walk constrained to return to its origin. 



2 Open Wilson loops (a toy example) 



We shall now proceed implementing the Wilson loop randomization program outlined in 
the introduction. To this end, we define the Brownian average of the Wilson loop as 



W Br (L) = I d£(s)5[F - 1] / dA; -trPexp 



L 2 /e 



ds^(s)A«t a 
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(14) 

The measure of the integration over paths £(s) can be straightforwardly understood: we 
are integrating over all functions £(s) subject to the constraint 5[C, 2 — 1] enforcing s to be 
the proper length parameter of the path. Note that we have replaced the total length of 
the path L by a formally divergent expression L 2 /e, where e is a cut-off associated with 
integrating over £(s). This is essential to maintain a meaningful continuum limit (as we 
shall see below explicitly). Indeed, random paths of a finite total length collapse to a point 
in the continuum limit. It is only paths of length L 2 /e that occupy a finite region of space 
(whose size is of order L) as e is taken to 0. This is immediately related to the fractal 
dimension of a random walk being equal to 2. 

The path integral over £ in ( I14p may appear complicated due to the non-linear constraint 
imposed by <5[£ 2 — 1]. This complexity is, however, spurious. The key idea for taking such 
path integrals over random walks is that the measure is completely factorized (<5[£ 2 — 1] 
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is a product of 5-functions at each point) and has the form of a factorized probability 
distribution in the celebrated central limit theorem of probability theory (with the limit of 
a large number of independent random variables being simply the continuum limit of the 
path integral). Then the result of the path integration (at least for observables depending 
on integrals of £ over finite intervals, rather than values of £ at a point) does not depend in 
the continuum limit on the precise form of the distribution of £, but only on the variance of 
this distribution. Hence, £[£ 2 — 1] can be replaced by a product of Gaussian distributions, 
which is straightforwardly treated in a path integral. 

There are many ways to implement the measure transformation we have just described. 
For example, [12] transforms the measure J d^(s)5[^ 2 — 1] using a Lagrange multiplier 
before invoking the central limit theorem. Here, we shall take a more straightforward 
route, manipulating the continuum limit of the measure in ( I14jl directly. (An alternative 
justification for all these formal manipulations with the path integral can be obtained, for 
example, by considering the continuum limit of random walks on a hypercubic lattice, which 
is an exactly solvable problem.) 

Let us first consider (for simplicity) the expectation value of a general observable de- 
pending on £ (integral of £), but not on the values of £ at a point (note that (TH|) is not 
quite of that form, but similar ideas apply, as we shall discuss below): 

(F) = J dZ(s)6[£ 2 - l]F(£( Sl ), • • • , £(**)). (15) 

(We shall assume Si < ■ ■ ■ < s/,.) Generically, one can expand F in a power series with 
terms of the form 

(^■)r---(^i)r- (is) 

Using f (s f ) = (f (s f ) - f (si_i)) + - £(si- 2 )) H h (f (si) - £(0)) (we assume that 

the random walk starts at £(0) = 0), it can also be re- written as a sum of terms of the form 




;, 



(£(**) - e(s fe -i)r • • • (£(si) - £(s )r = / dsn ■ ■• dsn . (17) 



Since the measure in (ITS]) is pointwise factorized with respect to £, and the factors in the 
last expression depend on £ (s) in non- intersecting ranges, the expectation value will likewise 
factorize (with an appropriate normalization of the path integral). It is therefore enough 
to consider the expectation value of just one factor: 

/ a > \ b 

2/ 



I(s,s') = / dv{s)5[v 2 {s) - 1] / dsv{s) 



where we have renamed £ to v . All the other observables of the type we are presently 
considering will be expressed as sums of products of such expectation values. 



5 



As was to be expected, (fT8j) is again of the form characteristic of the central limit 
theorem: a large (infinite in the continuum limit) number of uncorrelated variables with 
identical distributions, only appearing as an averaged combination in the observable to be 
studied. To study the continuum limit more explicitly, it is convenient to discretize the path 
integral in f lTSj) . splitting the interval (s, s') into (s' — s)/e points separated by distance e. 
One then obtains 

I £ (s,s') = ( d Vl l[5(v*-l)(eJ2^) b - (19) 

This integral will, of course, vanish for symmetry reasons unless b is even. Furthermore, 
when we expand {s^2vi) b , any terms containing odd powers of Vi will vanish after the 
integration. Hence, the structure of the expanded 6-power of the sum will be as follows: 
there will be ((s' — s)/e) b ^ 2 terms which are products of 6/2 values of vf at different points. 
For any terms with a higher power of Vi, one will obtain a smaller number of terms ( (s'—s) je 
to a power less than b/2) as more points have to coincide while picking v^s from the sum. 
Hence such terms will come with a greater power of e (than the terms containing only vf) 
and will not contribute to f lT9"j) in the e — > limit. Therefore, the leading small e behavior 
of (USD is 

I £ (s,s') = (e(s'-s)(v 2 )) b/ \ (20) 

where (v 2 ) = J dv v 2 5(v 2 — 1) = 1. (As already mentioned under ( TH|) . to obtain a continuum 
limit one needs to scale positions on the random walk as s ~ 1/e, which includes scaling 
the total length of the trajectory as L ~ l/e-) 

Note that ( 120]) is very insensitive to the shape of the distribution of v (as is typical of the 
situations in the spirit of the central limit theorem): any other distribution (different from 
8{v 2 - 1)) with the same (v 2 ) will give the same result for (I20p . and hence for expectation 
values of any quantities of the type ffT5j) . In particular, we can replace 5{y 2 — 1) by a Gaussian 
distribution exp(— v 2 ) (again, the normalization of the measure does not concern us here). 
In the continuum limit, ]^[exp(— vf) becomes exp (— J dsv 2 /e), which is in agreement with 
the continuum Gaussian measure obtained in [13J. The conclusion is that, for any quantity 
of the type (fT3|) . the random walk measure can be replaced as 



J di{s)5\i 2 -l}^ J d£(s)exv - J dsi 2 /e 



(21) 



without affecting the continuum limit (i.e., a random walk with a fixed speed becomes 
replaced by a walk with a Gaussian distribution of speeds, which is much more convenient 
to deal with in a path integral). 

Even though the randomized Wilson loop (THj) is not, strictly speaking, of the form (Tl5|) 
as the different integrations of £ are entangled with each other, a similar argument (though 
more cumbersome to typeset) leads to the conclusion that the integration measure in f)14p 
can be replaced as in (T2~T]) as well, without affecting the continuum limit. (The key idea 
is again that, after discretization, the only terms contributing to the continuum limit are 
products of vf at different points, and terms containing higher powers of v i are suppressed 
by powers of e. Hence, only the expectation value of v 2 contributes to the continuum limit, 
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and the (factorized) measure can be replaced by any other (factorized) measure with the 
same expectation value of v 2 . We use this freedom to replace the integration measure by 
the Gaussian one.) 

In this section, we are considering £(s) to be entirely unconstrained. In this case, the 
end point of the 'loop' does not have to coincide with the starting one. In a gauge theory, 
such a Wilson 'loop' would have to vanish by gauge invariance. However, in Eguchi-Kawai 
models this vanishing does not occur due to the breakdown of the R D symmetry. Hence, 
there is no straightforward relation between computations with closed and unconstrained 
Wilson loops (the situation may be better in more sophisticated reduced models, where 
open Wilson loops would have to vanish just as they do in gauge theory, and would thus 
give a zero contribution to the sum over all loops). Nevertheless, we shall pursue the 
computation of the unconstrained average to develop the relevant techniques, as this case 
enjoys a great algebraic simplicity. We shall follow up with a consideration of a sum over 
closed loops in the next section. 

Before we proceed evaluating the averaged Wilson loop (|14|) . we shall consider the 
leading large N large D behavior of 



J dA«e- s WA«\ ■ ■ ■ AX itr [t ai ■ ■ ■ t a2n ] . (22) 
Contracting all the A's according to ( {TBI , we obtain 

(g \ n 1 n 



{i(r)j(r)} 



Here, we have introduced pairings {i{r),j(r)} that for each pair number r give the two 
members of the pair i{r) and j(r), and the sum is over all such pairings. It is well-known 
from the usual considerations of the large N limit of gauge theories (see, e.g., [TJ [13]) that 
the leading large N contribution into tr [t ai ■ ■ - t a2n ] /N comes from planar contractions and 
it is equal to N n for any such planar contraction. 'Planar' in our context means that if all 
ta's are represented as dots on a horizontal straight line, and all the contractions <5 a - (r)aj(r) in 
( |23i) as arcs in the upper half plane connecting the corresponding points, such a diagram can 
be drawn without any intersections. The result can be represented pictorially as 'rainbow 
diagrams' or 'arches' of the type depicted in Fig. [TJ As only such diagrams are retained in 
the large N limit, one obtains 



planar{i(r)j(r)} 



gVN 



r=l 
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Figure 1: Example of a rainbow diagram. 



We now proceed evaluating the randomized Wilson loop. From (j!4p . §ZT§ and 

„ oo ^ l S ^ I 6 2n 

W Br (L) = j de(a)e-J' <fa ^ e X)(-ir J d Sl --- J ds 2n a M1 ... M2n J] £ 



S2n _! fc=1 

lMVeSr l_iXE\ [ / ^ 2n (25) 



dew c- w/e E l-^l I d Sl .- I ds 



n=0 



S2n-1 

2)i 



planar{«(r)j(r)} r=1 k=i 
The Gaussian integration over now merely amounts to contracting £ M 's with 

(e(si)e 2 (s 2 )) = £<^ 2 5(si - s 2 ). (26) 

In order to maximize the power of D in each term, ^'s should be contracted using exactly 
the same pairings {i{r),j(r)} that have been used for contracting A's in (122 p . With such 
a match, each S^kVi f rom the contraction (1261) will hit the corresponding £ Mi{r)Atj(r) in (I25p 
producing a factor of D. For any other contraction pattern, the number of powers of D 
will always be smaller, rendering such terms irrelevant in the large D limit. Bearing these 
observations in mind, we obtain 



oo / r-y L2 / £ l2 / £ 

W Br (L) = E I 2 I / ds l"" / rfs 2n E Yl S ( S i(r) ~ S Hr))i < 27 > 

™=° V / n o„.. , planar{i(r)j(r)} r=1 



where we have introduced the modified 't Hooft coupling A = g 2 ND. 

The rest of the computation becomes very simple now. The integration variables are 
arranged as s\ < s 2 < ■ ■ ■ < S2n- A 5-function 5(si( r ) — s_j( r )) enforces s^ r ) = Sj( r ) and hence, 
all Sj's with i between i(r) and j(r) should be equal to each other. Each such equality would 
make the integral vanish, except that the singularity in the 5-function can compensate for 
equality of two of the s^'s (but not more). That means that the integral in f[2TP will vanish 
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unless all the pairings are adjacent (j(r) = i(r) + 1). (This property can be derived more 
cleanly by considering a smeared (^-functions 5(s) = {1/2 A if \s\ < A, otherwise} with A 
that has to be taken to 0. One will then see that any integrals with non-adjacent pairings 
will be proportional to positive powers of A and thus vanish in the limit recovering the 
5-function proper.) 

There is evidently only one rainbow diagram with only adjacent pairings for each n. 
Furthermore, 

L 2 /e L 2 /e 

, n j f L 2 \ n 

dsf- I ds 2n Y[S(s 2r -i - s 2r ) = ^— j ( — J • (28) 



=i 





Hence, 

Wbt(L) = ^ — | -L\l'- ) =exp ( ) . (29) 

n=0 

So at the end of the day, all the dependences on the random path integration cut-off e have 
cancelled out (this simply reflects the fact that we have chosen the correct renormalization 
of the path length L = L 2 /e to maintain a finite spatial extent of the random walk), and 
we obtain a very attractive expression for the averaged Wilson loop (since L is the typical 
spatial extent of the random walk, L 2 is the typical area of its cross-section). 

As explained above, the derivations of this section are a toy example and do not bear 
an immediate relationship to the gauge theory (because the Wilson loops are not closed, 
and because of the spontaneous symmetry breakdown issues in the Eguchi-Kawai model). 
However, one may hope that the attractive mathematical structure we have displayed will 
survive in more realistic examples. Furthermore, in the next section, we shall repeat the 
present computation for the more realistic case of closed Wilson loops and derive (with 
considerably more effort) a somewhat more complicated formula resembling ( 1291) . 



3 Closed Wilson loops 

We shall now constrain the Wilson loop trajectory £(s) to start and end at 0. This results 
in a physically meaningful summation over gauge- invariant closed Wilson loops. 

Up to (1251) . the derivation is exactly the same as in the previous sectional Thereafter, 

1 Technically, for closed loops, £(s) are not completely independent at different points, as there is one 
global constraint f//° ds£ = 0. This constraint, however, does not upset the application of the central limit 
theorem. Indeed, one can make different £(s) completely uncorrelated by explicitly inserting S(J^° ds^) — 
Jdkexp (ikJ^dsiY With an unconstrained integration over £(s), one can consider again the expectation 

value of the type (IT5j) and include J dfcexp (ik j//° ds£j as a part of F. The derivations between (Tl5)) and 
(|2ip will then apply, the measure will become replaced by the Gaussian one, after which one restores 
S(J Q ds£) and removes it re-introducing the Dirichlet boundary conditions £(0) = £(Lq) = 0. The result 
is the original closed loop integral, but with the <5-function measure replaced by the Gaussian measure, as 
in (|21j) . Again, a variation of this argument will apply to the expectation value of the Wilson loop, which 
is not, strictly speaking, of the form (fl5]) . 
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to contract the £'s, one needs to use the Dirichlet Green function on a line: 



=eS"""D(s u s 2 ), (30) 
Dh,*)=ntoh,*)(l-2*p!»>). (31) 



(Once again, L = L 2 /e.) Or, for £: 



(i^{8i)i^{82)) = eS^V{ 8l ,8 a ), (32) 
V(s 1 , s 2 ) = d sl d S2 D(s 1 , s 2 ) =5(s!- s 2 ) - — . (33) 



Then, 



W^r(^) = Xl(- £ V^) / / ^ 2n ^ II^(^(r),^M)- (34) 

«=0 V / S2n _ 1 planar{i(r)j(r)} r=l 

This expression only differs from f l27|) by the second (constant) term on the right-hand 
side of (13"3"|) . Yet, such a minor difference turns out to complicate the remaining derivation 
considerably. We shall first re-scale s so that it runs from to 1: 

w Br (L) = jr | -l 2 ^| J Jd Sl ■ ■ ■ J ds 2n n (%(»•) - - x ) • 

n=0 \ 'a S2 „_! planar{«(r)i(r)} r=l 

(35) 

One sees that the random walk integration cut-off £ disappears as it should. Expanding the 
product in the above formula, one will have a collection of terms with different numbers 
of ^-functions. As explained under ( 1271) . these integrals will vanish, unless all the pairings 
appearing in the ^-function are adjacent (j(r) = i(r) + 1). If that is so, a general term of 
this type will have k ^-functions whose argument is a difference of adjacent Sj's (denoted 
s p (i) and s p (/)_|_i), and it can be evaluated as 

i i , 

Jdsf- J ds 2n f[ 5{s p{l) - a p(Q+1 ) = 2k , 2 *_ k y ■ (36) 

s 2n -l ' =1 

We now introduce i? mn , the number of rainbow diagrams (Fig. [I]) with m adjacent and 
n non- adjacent pairings, and consider a rainbow diagram with a total of n points, m of 
which are adjacently-paired. For all non-adjacently-paired points, the 5-function term in 
( 1351) cannot contribute. For the adjacently-paired points, one can have contributions from 
both ^-function and constant term. Furthermore, for k ^-functions, there are ways to 
choose their positions on the m arches connecting adjacent points, all resulting in the same 
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value of f )36|) . Hence, f )35|) can be re- written using f )36|) as 



00 / / \ \ n m / 1 \n-kr<k 

«wl)=E(-W£) E^E^bt 

n=0 \ / m=0 fc=0 V 7 




What remains now is to evaluate the rainbow numbers R mn and to resum ([3 
The rainbow numbers satisfy a recursion relation 

m+ 1 n 



(37) 



+ R klRm+l-k,n-l (38) 



fc=l i=0 

with the boundary conditions 

i? n0 = 1, Rm = 1 (any n), R 0n = (n > 0). (39) 

This recursion relation can be easily understood: consider a rainbow with m + 1 adjacent 
and n + 1 non-adjacent pairings. If the first point is connected to the second one (i.e., 
adjacently paired), the remaining diagram is a rainbow with m adjacent and n + 1 non- 
adjacent pairings. This gives the first term. Alternatively, the first point can be connected 
to another point with an even number. Then, the arch coming out of the first point will 
split the diagram into two smaller rainbow diagrams (under this arch and outside this arch) 
with a total number of m + 1 adjacent and n non- adjacent pairings. This gives the second 
term. 

The recursion relation (1381) is closely reminiscent of the one for Catalan numbers. This 
is not surprising, since the total number of rainbows R n = Yl m Rm,n-m is precisely given 
by Catalan numbers (see p. 222-223 of [1]). To solve the recursion relation, one introduces 
a generating function 

00 

R(x,y)= R mnx m y n , (40) 

m,n=0 

which allows to rewrite f)38p and f)39p as 

1 ( R 1 ) = i( R 1 ygJ. (41 ) 



xy \ 1 — x J y \ 1 — x J x 
This quadratic equation is solved by 
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Keeping (HOI) in mind, we must choose a solution non-singular at the origin, which is the 
one with a minus sign. This non-singular solution can also be re-written as 

R = 2 . (43) 

1 — x + y + — x + y) 2 — Ay 



As a cross-check, the generating function for the Catalan numbers (1 — a/1 — 4x)/2 should 
result from setting x = y, which is indeed the case. The first few terms of the power series 
for R are 



R = 1 + x + x + xy + x + 3x y + xy + x + 6x y + 6x y + xy 

(44) 

+x 5 + 10x 4 y + 20x 3 y 2 + 10x 2 y 3 + xy 4 + ■ ■ ■ , 

and they can be verified to match the numbers of low order rainbow diagrams. 

A brute force way to proceed with the computation would now be to extract the rainbow 
numbers by expanding R and then try to resum (1351) . However, one can do wiser than that 
and attempt to convert the (already resummed) expression for R into an expression for the 
Wilson loop. To this end, we observe that 

OO 771 , 

*(*(l-f),y)= £ R mn x m y n C k m (-£) (45) 

m,n=0 k=0 

and 



R (f i 1 ~ I) • JO = T, «-*v E ct (-5) ^ 2 "" + " ,+t . («) 

Furthermore, using 



m,n=0 k=0 



2m J z n+1 n\ 

c—io 

(where c is a positive number), we obtain 



c+ioo 

^7 I dz4- = ^ (47) 



e+ioo ^ 

— / ds-flf-^fl—V-^ ^"yV w l" 1 ^ (48) 
2vri ./ z \z 2 \ 2J , z 2 ) ^ y ^ 2 k (2(n + m) - k) K J 



_ m,n=0 k=0 

Here, we immediately recognize a structure matching fl35l) . Hence, with £ 2 = L 2 a/A/2, 

c+ioo 



(49) 

C+IOO 

If. 4e 2 
/ 

2vri 



2z + r + Ay(2^ + ^) -16£ 2 
12 



Because of the fast decay of the exponential, the contour can actually be closed at z = — oo 
and then pulled over the cut of the square root in the integrand. We obtain 



c-Moo 2£ 

„, , r , exp(-£ 2 /2) f , e z exp {-£ 2 /2) f, x r— 

Wbt(L) = . / dz . = ' ' / dx eV4£ 2 - x 

™ J z + Vz 2 - U 2 2tt£ 2 J 

c—iao —2£ 



(50) 

2 V 2 / 7T 



expl-^/rM - / da e 2L ^ lfi — sin 2 a. 



o 



The leading behavior here is again in the form of an area law: 



4 Conclusions 

Considering randomized Wilson loops (Wilson loops averaged over an ensemble of random 
contours) in the context of Eguchi-Kawai reduced model results in an algebraically tractable 
dynamically non-trivial large D limit (D being the number of matrices in the Eguchi- 
Kawai model, corresponding to the number of dimensions in gauge theory). The scale 
dependence of the expectation values of such randomized Wilson loops turns out to be 
strongly reminiscent of the area law. 

It would be interesting to generalize the present considerations for randomized Wilson 
loops to more realistic reduced models [2], [7] . This should be possible once the large D limit 
has been worked out for those theories. The area law seems to emerge quite naturally from 
the algebraic structures appearing in the randomized Wilson loop computations (a series 
in powers of L 2 , i. e., the area, with alternating signs and factorial denominators). It will 
not be surprising if analogous computations for more realistic cases display an area law as 
well. 

It is also interesting to speculate how much the behavior of randomized Wilson loops 
constrains the behavior of individual Wilson loops. The area law decrease for large Wilson 
loops is very fast, and having individual Wilson loops decreasing slower than the area 
law would seem in conflict with the fast decrease of the randomized average. Should it 
be possible to show that expectation values of individual Wilson loops are positive, the 
fast decrease of the randomized average will become even more constraining (if a sum of 
positive numbers is small, it constrains the magnitude of individual entries in the sum). 
A better understanding of the relations between randomized and individual Wilson loop 
averages could also help to establish the randomized Wilson loop analyzed here as an order 
parameter for the confining phase of the gauge theory (cf. the preceding remarks on the 
area law behavior). These issues certainly deserve further investigation. 

Since the matrix integral ([I]) appears as a small volume limit of compactified gauge the- 
ory (irrespectively of the Eguchi-Kawai equivalence), our present derivations formally apply 
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to the study of that limit as well. The caveat is that the matrix model Wilson loop (jlj) does 
not appear to have special significance, unless it can be related to the gauge theory Wilson 
loop via the Eguchi-Kawai equivalence. Some applications of the high temperature/small 
volume limit of gauge theories can be found in [T4] . 

It is instructive to examine our results in the light of numerical findings reported in [15J. 
In that publication, area law for Wilson loop expectation values has been seen via numerical 
simulations for sufficiently large rectangular Wilson loops in bosonic and supersymmetric 
matrix models similar to (JTJ). This observation was perceived as an evidence that the 
Eguchi-Kawai equivalence holds for sufficiently large spatial scales (since the corresponding 
gauge theories are confining and should indeed exhibit the area law). The validity of 
the Eguchi-Kawai equivalence for sufficiently large scales has been further heuristically 
linked to another numerical result described in [15] : the smallness of Polyakov lines of 
sufficiently large circumference. Since vanishing of Polyakov lines would imply a lack of 
spontaneous breaking of the U(N) center symmetry and hence exact validity of the Eguchi- 
Kawai equivalence, the observed smallness of Polyakov lines was seen as an additional 
indication of approximate validity of the Eguchi-Kawai equivalence for sufficiently large 
spatial scales (physically, above the characteristic size of confined one-particle states). All 
these observations mesh quite well with the derivations presented here: if the heuristics 
of [15] is taken one step further, one may hypothesize that the Eguchi-Kawai equivalence 
applies to randomized Wilson loops of sufficiently large characteristic size (L). Then the 
present results will have direct implications for gauge theory irrespectively of spontaneous 
breaking of the U(N) center symmetry that undermines exactness of the Eguchi-Kawai 
equivalence. Of course, such claims are far from high standards of rigor, as are those of 
[15J. Whatever the case, the derivations presented here can be seen as a first step towards 
analytic understanding of the numerical observations reported in [IS] . 
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